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Abstract 

Background: The signals that the central nervous system (CNS) produces and sends to the muscles to effect 
movement are not entirely understood. Muscle synergy theory suggests that the central nervous system produces a 
small number of signals that pass through a network that distributes combinations of these signals to the muscles. 
Though these synergies are rather stable over time, some variability is present. 

Methods: Here, we investigated the variability of muscle synergy and defined a synergy stability index (SSI) to 
quantify it. We measured the activity of muscles responsible for maintaining lateral balance in humans standing on a 
platform that was subjected to lateral disturbance from the platform. We then calculated muscle synergies attributed 
to postural reflex and automatic response by using non-negative matrix factorization (NMF). Finally, from the 
calculated muscle synergies, we obtained SSI. 

Results: We observed a positive proportional relation between balance performance and SSI. Participants who were 
adept at maintaining balance were found to have invariant muscle synergies, and non-adept participants showed 
variable muscle synergies. 

Conclusions: These results suggest that SSI can be used to quantitatively evaluate balance maintenance ability. 



Background 

The term muscle synergy, also called motor synergy, neu- 
romuscular synergy, or muscle mode, has been used in 
the literature extensively over the last decade. A muscle 
synergy is the activation of a group of muscles to con- 
tribute to a particular movement [1], thus reducing the 
dimensionality of muscle control A single muscle can be 
part of multiple muscle synergies, and a single synergy 
can activate various muscles. Different ways of group- 
ing muscles into synergies can be found in the literature 
[2]. Some researchers define synergies as muscle acti- 
vation of a set of muscles contributing to a particular 
movement where each muscle contributes to only one 
synergy [3]. However, the existence of synergies is con- 
troversial, with studies providing evidence for and against 
the existence of muscle synergies [2]. Muscle synergies 
are calculated from electromyography (EMG) data, and 
have been described in frogs [4], cats [5], and humans 
[3,6-8]. Studies have focused on voluntary movements 
such as walking [9], cycling [10], and upper-limb reaching 
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[6]. Others have investigated reflexes in frogs by cuta- 
neous mechanical stimulation [11] and in cats through 
displacement of a supporting surface [5], as well as auto- 
matic postural responses in humans, where upright stance 
synergies were identified for different directions [7] . 

Muscle synergies can be identified by factor analyses 
[12] such as non-negative matrix factorization (NMF) 
[13,14] and principal component analysis (PCA) [15]. 

In this study, we calculated synergies from EMG data 
and used an NMF method based on the Lee-Seung algo- 
rithm [14], which is a basic and fast NMF algorithm with 
multiplicative updating rules. 

We chose balance maintenance for our task to observe 
reflex and automatic responses of the central nervous sys- 
tem (CNS). We based our study on our previous work 
on human balance [16]. In a society where humans are 
living longer than ever before, deterioration of ability to 
maintain balance poses a serious problem for elderly per- 
sons [17]. The mechanism of balance maintenance must 
be clarified in order to find a solution to this growing 
problem. 

Our hypothesis about the control architecture of the 
CNS is based on the bow-tie structure proposed in [18]. 
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We regard this hypothesis as an assumption for the exis- 
tence of synergies. Only a portion of the information 
about the environment outside the human body enters 
the sensory system. An even smaller portion of this infor- 
mation is relevant to a particular task, such as balance 
maintenance. The CNS processes this relatively small 
amount of information and sends signals to the muscles 
that, in turn, move the highly redundant musculo-skeletal 
structure. The bow-tie structure represents a model of 
reducing degrees of freedom, since it processes high- 
dimensional input and produces high-dimensional output 
by a rather low-dimensional processing unit. Dealing with 
low-dimensional data is very time-efficient and allows for 
processing multiple tasks at the same time. 

We cannot measure the relevant synergy recruitment 
signals directly at source, that is, in the brain or other parts 
of the CNS. However, we can measure the outcome which 
is the muscle activation. From these data, we can estimate 
what signals the CNS must have sent. We call these signals 
synergies. 

Since all muscles move the same body, and there are 
multiple muscles moving the same joint, most muscles 
work in groups. The muscle activation signals are related 
to each other and thus they can be represented by a 
smaller number of signals. From the EMG signals, we cal- 
culate these representative signals (synergy recruitment 
signals) and the weights (muscle synergies) of the net- 
work that distributes these signals to the muscles. The 
number of synergy recruitment signals and the number of 
synergies are therefore much smaller than the number of 
muscles used. 

Here, we propose an index of synergy stability to quan- 
tify synergy variability. Our findings suggest that synergy 
stability can be a measure of task performance quality 
and can be quantified by using the metrics introduced 
hereafter. 

Methods 

Definition of muscle synergy 

We assume that each EMG signal Mj (t) from a muscle 
that produces a reaction to a disturbance of body posture 
can be described as 

k 

Mj(t) = J2wijQ(t)> j = l, ••• ,m (1) 

i=l 

where m denotes the number of measured muscles, Q (t) 
is the signal produced by the CNS to control the /th syn- 
ergy, Wij is a time-independent weight for muscle j in 
synergy /, and k is the number of synergies. We define 
Mj (t) for the interval from 75 ms to 150 ms after dis- 
turbance onset. We chose this interval because reflex 
mechanisms and automatic responses are active within 
this time window. 



Lateral disturbance experiments 

All experiments were approved by the RIKEN Ethics 
Committee. 

Ten healthy persons (all men; age, 22-50 years; weight, 
55-95 kg; height, 164-185 cm) participated in the exper- 
iment (Table 1). Participants stood on a platform with 
arms akimbo (on their hips) and feet parallel at about 30 
mm apart, as shown in Figure 1. The arms akimbo posi- 
tion was used to cancel the effect of balancing using the 
arms. It is also closer to the natural hanging-down posi- 
tion than crossing arms on the chest, which made many 
participants feel uncomfortable. Crossed arms would also 
have obstructed optical markers placed on the chest for 
motion capture. Distance between the feet was set slightly 
shorter than shoulder width because participants did not 
show balance difficulties at greater distances. To effec- 
tively disturb the balance of a person with feet set at 
shoulder-width distance, the disturbance would have to 
be much stronger, which was technically infeasible and 
potentially dangerous. 

During each trial the platform moved laterally as shown 
in Additional file 1. The timing and direction of move- 
ment were random and not known to the participants 
in advance. Movements in only the left direction were 
used for evaluation. During the movement of the plat- 
form to the left, the human body inclines to the right, 
putting body weight almost entirely on the right leg. 
The participants were instructed to stand still, keep their 
feet on the ground and avoid stepping and body move- 
ments other than lateral hip flexion/extension and ankle 
inversion/eversion. EMG signals of six muscles of the 
right side from 75 ms [19] to 125 ms after movement 
onset were recorded by using wireless active electrodes 
(BTS Bioengineering, Milan). The six muscles chosen 
for EMG recording were relatively easy-to-measure sur- 
face muscles in charge of lateral posture correction: 
muscles for ankle joint inversion, the flexor hallucius 
longus (FHL) and tibialis anterior (TA); for hip joint 



Table 1 Participant data 



Participant 


Age 


Height [cm] 


Weight [kg] 


sbjl 


33 


183 


72 


sbj2 


40 


172 


68 


sbj3 


50 


164 


63 


sbj4 


28 


185 


95 


sbj5 


45 


170 


59 


sbj6 


35 


170 


60 


sbj7 


33 


172 


72 


sbj8 


30 


178 


73 


sbj9 


49 


183 


79 


sbjl 0 


22 


165 


55 
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Figure 1 Participant standing on the platform. 



abduction, tensor faciae latae (TFL) and gluteus medius 
(GM); and for trunk lateral flexion, rectus abdominis 
(RA) and erector spinae (ES). Because RA produces a 
noisy signal when there is a fat layer over the abdomi- 
nal muscles, only participants with relatively low fat were 
chosen. 

The following measures were implemented to create 
equal conditions for all participants and to keep them as 
motivated as possible. 

We introduced a score system (see the subsection on 
assessment of inter-individual differences) to help partic- 
ipants maintain their concentration and motivation. The 
system provided a platform score (PS) that was displayed 
on a screen to give the participants feedback on their 
performance. 

The participants stood on thin force sensors used to 
monitor lateral balance. The experiment was executed 
after confirming that the participants were balanced, that 
is, center of pressure was at the midpoint between the feet. 
This was done to prevent the participants from trying to 
predict the direction of the platform movement. 

Muscle activation was monitored, and the experiment 
was carried out only when EMG activity was less than 
30% of maximum, indicating a relaxed posture. This was 
done to prevent the participants from trying to predict the 
timing of the platform movement. 

Two different disturbance strengths were used. The 
platform moved 160 mm within 0.66 s for the weak dis- 
turbance and 180 mm within 0.75 s for the strong dis- 
turbance. The parameters (displacement and duration) 
of the weak disturbance were tuned so as to produce a 
muscle response strong enough to obtain an acceptable 
signal to noise ratio but weak enough to not evoke a 
stepping response. The strong disturbance was used to 
determine differences in behavior between adept and non- 
adept participants for balancing. Adept participants were 
able to maintain their balance, while non-adept partici- 
pants responded to the strong disturbance with a step or 
by lifting a leg. 

Calculation of muscle synergies from EMG data 

EMG was sampled at 1 kHz, high-pass filtered with a 
30 Hz cutoff, root-mean square rectified, and smoothed 
using moving average with sine function shaped window 
of length 10 samples. 

The time window used for synergy extraction from 
EMG signals varies between studies. Many researchers 
investigating automatic postural responses base their 
methods on [20], where averages were taken from three 
time bins, two 75 ms long followed by one 350 ms long, 
starting 75ms after onset of a disturbance. Similarly in [7], 
the authors used three 75-ms time bins after disturbance 
onset. A finer time division was also applied for voluntary 
movements in [8]. 
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Muscle activity signals Mj (t) (Eq. 1) constituted the 
experimental data and were arranged to form an m x n 
matrix M, where m denotes the number of measured mus- 
cles and n the number of samples. EMG data M were then 
factorized by NMF to obtain the synergy matrix W and 
the synergy recruitment matrix C such that 



(2) 



where k denotes the number of synergies and E the resid- 
uals. EMG data from each trial were used as input to the 
NMF algorithm. We selected the Lee-Seung algorithm 
[14] for our NMF, as other algorithms tended to be unsta- 
ble. The initialization matrix of the NMF algorithm was 
set randomly, and the calculation was run 10 times for 
each data set. No significant changes in the results were 
observed. 

Each of the 10 participants performed total of 12 left- 
direction and 12 right-direction trials in random order. 
Four left-direction trials per person were selected for fur- 
ther synergy calculation. Those were trials where the par- 
ticipants did not expect the disturbance and were relaxed. 
Specifically, they showed balanced posture, equal values 
for the foot sensors, and low EMG activity shortly before 
the onset of the disturbance. On the basis of observations 
and questioning of the participants, we decided to treat 
the first trials as a warm up and excluded them from eval- 
uation. Data from one participant who reported fear and 
fatigue on the day of the experiment and performed many 
unnecessary recovery movements were discarded as out- 
liers. Data from the remaining 9 participants were used for 
evaluation. 

Number of synergies 

We calculated the variability accounted for (VAF) to 
determine to what extent the original EMG data could 
be reconstructed from the factorized data. The VAF is 
defined as 



VAF = 1 - 



\\E\\ 2 F 
\\M\\ 2 F 



2 ' 



(3) 



where ||-||p denotes the Frobenius norm of a matrix. VAF 
calculations were carried out for each number of syner- 
gies 1 to m. When VAF was higher than the 90% threshold, 
then the number of synergies was deemed sufficient to 
regard the factorization results W and C as representative 
of the data setM [21,22]. 

Assessment of inter-individual differences 

To assess inter-individual differences, we introduced two 
quantitative indexes of balance ability, synergy stabil- 
ity index (SSI), based on synergy calculations and plat- 
form score (PS), based on observations of human balance 
recovery. 



Synergy stability index 

To obtain an SSI for each participant, synergies were cal- 
culated for each weak disturbance trial, and their variabil- 
ity was evaluated by calculating the correlation between 
the synergies of that participant. 

The rows of the matrix W are the synergy vectors, and 
we define wj (j = 1, . . . ,p) as the synergy vector of a par- 
ticipant for the yth trial for a particular synergy. If wj are 
well aligned, or more precisely, if w\ = Wi — ... = 
w p = w after normalization, then we can say that w 
is the synergy vector of this participant for the experi- 
ment. In reality, synergy vectors are not identical from 
trial to trial, and we evaluate the degree of alignment 
among synergy vectors of different trials, representing this 
by the Pearsons correlation coefficient r (Eq. 5) for each 
of the synergy vector pairs. Results were averaged over all 
p (p — 1) /2 pairs, and then the averaged pairwise corre- 
lations were further averaged with respect to synergies to 
yield SSI. More precisely, SSI is represented as 



1 k 

ssr- £ 



i=l 



p(p- 1) 



(4) 



where p is the number of trials, k is the number of syner- 
gies, and and vJq are the /th normalized synergy vec- 
tors of the /th and #th trials, respectively. Synergy vectors 
were normalized such that the sum of all elements squared 
equaled 1. We sorted the resulting synergies to determine 
the highest SSI. The Pearsons correlation coefficient r is 
defined as 



r(x,y) = 



E£i(*y-*)(yy-5>) 



W-SxSy 



(5) 



where x and y are two vectors to compare; x and y are 
mean values, while S x and S y are standard deviations. 

Platform score 

Balance ability was also assessed by observing participants 
behavior during balance recovery from the strong distur- 
bance and quantitatively expressed by PS. All trials were 
used for calculation so as to evaluate the overall ability of a 
participant and not a single strategy as in SSI calculation. 

The maximum score was achieved when participant 
stood on both feet and kept both hands in the akimbo 
position at all times during the disturbance. Scores were 
penalized for the following behaviours. Removing one 
hand from the akimbo position was penalized one point, 
and removing both hands from the akimbo position was 
penalized two points. Lifting one heel off the surface was 
penalized one point, completely lifting one foot off the 
surface was penalized two points, and stepping aside with 
one foot was penalized three points. Totally losing balance 
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Figure 2 VAF.Two examples of VAF. SSI for number of synergies k = 1 ... 6 are shown. Participant sbjl (left) has one synergy since his SSI for k > 1 
is greater than the 90% threshold (horizontal line), and participant sbj9 (right) has two synergies since his SSI for k > 2 is greater than the threshold. 
Error bars denote standard deviation. 





FHL TA TFL GM RA ES 0 20 40 60 

muscles time [ms] 



Figure 3 Synergy. Two sets of synergies W (bar graphs) and synergy recruitments C (time plots) for participant sbj6. a: weak disturbance; b: strong 
disturbance. 
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and stepping aside with both feet cost the participant all 
six points. 

Results and discussion 

VAF calculation (Figure 2) revealed that the number of 
synergies k needed to reconstruct the original muscle 
activity data of all measured muscles was one or two, 
depending on the participant. Thus, we extracted an equal 
number of synergies (two) across participants to facili- 
tate comparison. This low number of synergies can be 
explained by the fact that the human body in the lateral 
direction can be modeled with one degree of freedom for 
weak disturbances and two degrees of freedom for strong 
disturbances. For weak disturbances, the ankle joint is 
sufficient to maintain balance (ankle strategy); for strong 
disturbances, the ankle together with the hip joint play the 
main roles in balance maintenance (hip strategy) [23] . 

Figure 3a shows an example set (sbj6) of synergy W and 
synergy recruitment C. Each of the two synergy vectors W 
can be found in a separate bar graph on the left side. Each 
of the two synergy recruitment signals C is presented in a 
separate graph on the right side. 

Mean EMG values and the standard deviations of all six 
measured muscles (Figure 4) revealed that muscle acti- 
vation can vary from trial to trial. The synergy similarity 
among trials is expressed by the SSI. Its values, for each 
participant, can be seen in Figure 5b. However, despite 
inter-trial differences among EMG signals, inter-trial syn- 
ergies revealed to be similar as expressed by SSI values 



(Figure 5b). PS mean and standard deviation (Figure 5a) 
were calculated from approximately 12 platform move- 
ments for each participant and provided a quantification 
of balance performance. PS was not found to be depen- 
dent on height or weight and to be only slightly dependent 
on age (Figure 6). 

Plotting PS as a function of SSI for each participant 
(Figure 5c) and applying linear least squares regression 
revealed a strong correlation (R = 0.82) for the nine par- 
ticipants. This indicates that participants adept at main- 
taining balance (high PS) used the same synergy (high 
SSI), while non-adept participants used different syner- 
gies (low SSI) across trials. In other words, it can be said 
that persons who have high PS also have high SSI, and 
vice versa, as can be seen in Figure 5c and by comparing 
Figure 5a with Figure 5b. Therefore, participants who were 
able to minimize synergy variability were the most adept 
at maintaining balance. 

PS was derived from strong disturbance trials to allow 
observation of differences in balance ability between par- 
ticipants. The strength of the disturbance was chosen such 
that some participants were not able to maintain their bal- 
ance without waving their arms or stepping aside. SSI, in 
contrast, was derived from weak disturbance experiments. 
At this strength of movement, all participants were able 
to maintain balance without arm movements or stepping. 
One participant who was adept at balancing for both weak 
and strong disturbances showed similar synergies for the 
two disturbances (Figure 3a and b, left). Differences in 
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Platform Score PS 




CO 
Q_ 



sbj2 sbj5 sbj9 sbj10 sbj4 sbj6 sbj7 sbj8 sbjl 
subject 

Synergy Stability Index SSI 




sbj2 sbj5 sbj9 sbj10 sbj4 sbj6 sbj7 sbj8 sbjl 
subject 

regression: SSI - PS, R=0.823 




0.4 0.5 0.6 
SSI [1] 

Figure 5 PS and SSI. a: PS of nine participants from the strong 
disturbance; b: SSI of nine participants from the weak disturbance; 
c: Relation between SSI and PS, with linear least squares regression 
(solid line). 



Platform Score (PS) - age, R=-0.688 




4.5 
4 

3.5 
3 

2.5 

CO 

Q_ 2 
1.5 
1 

0.5 

0 



4.5 
4 
3.5 

3 

— 2.5 
(I) 

Q_ 2 
1.5 
1 

0.5 



5 30 35 40 i 

age [a] 

Platform Score (PS) - height, R=0.228 




165 170 175 180 185 

height [cm] 



Platform Score (PS) - weight, R=0.167 







o 








o c 


> 












o 




o 








o 

I 1 





50 60 70 80 90 

weight [kg] 

Figure 6 PS as function of age, weight, and height. PS as a 

function of age, weight, and height, with linear least squares 
regression (solid line). 



Wojtara et al. Journal ofNeuroEngineering and Rehabilitation 2014, 1 1:1 29 
http://www.jneuroengrehab.eom/content/1 1/1/129 



Page 8 of 9 



muscle activity (Figure 3a and b, right) can be attributed 
to the synergy recruitment C. 

The results suggest that the CNS is not able to produce 
a constant synergy in persons who are not adept at main- 
taining balance. Adept persons, on the other hand, will 
show less variations between disturbances. We assume 
that during balance training the CNS tries to adjust its 
synergy within a certain range to find the best configura- 
tion. In non-adept persons, this synergy is adjusted in a 
wider range than in persons who are adept at maintaining 
balance. 

Our study was not without limitations. Our balance 
scoring system was not a clinically established index of 
balance. However, we were forced to use a self-designed 
scoring system because established indexes didn't fit our 
needs. The Berg Balance Scale, for example, is a well 
established index that neglects dynamic balance and was 
designed for persons with balance disabilities while we 
dealt with healthy persons and needed an index to exam- 
ine their dynamic balance. Furthermore, EMG data for 
the ES muscle tended to have noise due to sweat on the 
lower back as a result of physical exertion during the 
experiment. 

Conclusions 

We applied a platform disruption method to investigate 
the neural mechanisms of involuntary motor control. 
Muscle synergies with distinct space-time separation were 
found, suggesting a possible method for reduction of mus- 
cle redundancy in motor control. We measured muscle 
activity of reflex and automatic posture control, then com- 
puted the associated muscle synergies. Analysis revealed 
that some participants showed relatively constant syn- 
ergy across trials, while others showed different synergies. 
We demonstrated that participants with constant syn- 
ergy were more adept at maintaining balance after lateral 
perturbation than those with different synergies among 
trials. This suggests the potential for the use of synergy 
constancy represented by SSI as a measure of ability to 
maintain balance or neural control quality. 

Consent 

Written informed consent was obtained from the patient 
for the publication of this report and any accompanying 
images. 
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